Pion dominance in I^p SUSY induced neutrinoless double beta decay 
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At the quark level there are basically two types of contributions of R-parity violating SUSY 
(-^p SUSY) to neutrinoless double beta decay: the short-range contribution involving only heavy 
virtual superpartners and the long-range one with the virtual squark and neutrino. Hadronization 
of the effective operators, corresponding to these two types of contributions, may in general involve 
virtual pions in addition to close on-mass-shell nucleons. From the previous studies it is known that 
the short-range contribution is dominated by the pion exchange. In the present paper we show that 
this is also true for the long-range I^p SUSY contribution. Therefore, we conclude that the I^p SUSY 
contributes to the neutrinoless double beta decay dominantly via charged pion exchange between 
the decaying nucleons. 
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I. INTRODUCTION 

The nuclear neutrinoless double beta decay {01/(3(3) is 
a process known for 70 years, which has been searched 
for, but not yet seen. The Oi^/?/3-decay plays a special 
role among exotic processes for the following two rea- 
sons. First, it is intimately related to the nature of the 
neutrino since it is able to probe whether the neutrino 
is a Dirac or a Majorana particle. Second, the modern 
0^/?/3-experiments have reached unprecedented sensitiv- 
ity. The most stringent lower bound on the half-life of 
the Qi/(3P-decay was measured for ^^Ge in the Heidelberg- 
Moscow experiment [2:] : 



Ti/'2{ Ge) > 1.9 X lO^Vrs 



(1) 



In the near future this limit is expected to be improved 
in the GERDA experiment |3| by 1-2 orders of magni- 
tude. For the Oj//3/3-decays of i°°Mo and "°Te in the 
two running experiments, NEMOS [i| and CUORICINO 
0, the following sensitivities have been achieved: 



T°;2('°°Mo) > 5.8 10^3 years 
T^J^C^^Te) > 3.0 10^4 years. 



(2) 



The extraordinary sensitivity of 0^/3/?-experiments 
makes them also a unique laboratory tool to probe 
physics beyond the standard model (SM) with possible 
lepton number violation (LNV) underlying. Of partic- 
ular interest in this context are supersymmetric models 
with R-parity violation {^p SUSY) containing the LNV 
interactions necessary to generate the Majorana mass 
for neutrinos and to induce the Oj^/3/3-decay. There is 
a wealth of literature on the R„ SUSY mechanisms of 
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the 0i'(3(3-dec<iy d S i, i, [O, E El, Q [ll where 
the corresponding quark-level LNV interactions as well 
as hadronic and nuclear structure aspects relevant for 
these mechanisms are studied. At the quark-level there 
are basically two types of J^p SUSY mechanisms: the 
short-range mechanism with the exchange of heavy su- 
perpartners [1, 0, S, M, [13, [m and the long-range mech- 
anism involving both the exch ang e of heavy squarks and 
the light neutrino [l^, [H, 0, Il5l | which we call squark- 
neutrino mechanism. For the latter case, due to the chi- 
ral structure of the ^p SUSY interactions, the amplitude 
of Oi^/3/9-decay does not vanish in the limit of zero neu- 
trino mass in contrast to from the ordinary Majorana 
neutrino exchange mechanism proportional to the light 
neutrino mass. Instead, the squark-neutrino mechanism 
is roughly proportional to the momentum of the virtual 
neutrino, which is of the order of the Fermi momentum of 
the nucleons inside the nucleus with pp ~ lOOMeV. This 
is a manifestation of the fact that the LNV necessary for 
Oj^/3/3-decay is supplied by the I^p SUSY interactions in- 
stead of the Majorana neutrino mass term and, therefore, 
this mechanism is not suppressed by the small neutrino 
mass. 

In the calculation of the Oi^/3/3-decay amplitude one 
has to hadronize the quark-level LNV operators repre- 
senting them in terms of the interpolating nucleon fields. 
As is known (9|, this can be done in two ways: the quark 
fields can be directly imbedded into the interpolating two 
nucleon fields (the two nucleon mode) or via an interme- 
diate step when one quark-antiquark pair is associated 
with the charged pion field coupled to the nucleons (pion 
mode). In Refs. %^ it was shown that the pion mode 
of hadronization dominates over the two nucleon one for 
the case of the short-range mechanism. In the present pa- 
per we demonstrate that the same conclusion is valid for 
the long-range squark-neutrino mechanism. Therefore, in 
all cases the dominant contribution of the -^pSUSY to the 
Oj^/?/3-decay is realized via the pion mode of hadroniza- 
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tion. For the calculation of the nuclear matrix elements 
(NMEs) of the corresponding transition operators we ap- 
ply the proton-neutron Quasiparticle Random Phase Ap- 
proximation (pn-QRPA) [l^. 

The paper is organized as follows. In Sect. II we in- 
troduce the quark-level process related to the squark- 
neutrino mechanism of the Oiyf3(3. In Sect. Ill we discuss 
the two nucleon and the pion modes of hadronization of 
the effective quark-level Lagrangian. The NMEs of 

the corresponding transition operators are presented in 
Sect. IV. For the three nuclei of experimental interest the 
values of these NMEs are calculated in the QRPA. They 
are discussed in Sect. V. We summarize and discuss our 
results in Section VI. 



II. SQUARK-NEUTRINO MECHANISM OF 

Oi//3/3-DECAY 

Lepton number violation in the fermions sector is gen- 
erated within SUSY models by the following terms of 
the superpotential: 

Wijt^ = ^KjkULjEk + X'.jkLzQjDk + n'UH2, (3) 

where L and Q are the left-handed lepton and quark su- 
perfield SU{2) doublets, while E, U and D denote the 
right-handed lepton, up-quark and down-quark SU (2) 
singlets, respectively. 

The mixing between the scalar superpartners Ql.r. of 
the left and right-handed quarks qL,R plays the crucial 
role for the case of the long-range SUSY squark- 
neutrino (SQN) mechanism. This effect occurs due to 
the non-diagonality of the squark mass matrix. For the 
down squarks of each generation it takes the form (e.g. 

0) 



a 



9 

m r 



-md{Ad 



- QA2Dz 
/i tan/3) 



-■md{Ad 



/i tan/9) 



(4) 



Here, d = d, s,b and d are their superpartners. Dz = 
M|cos2^ where tan^ = {H^)/{H^) is the ratio of vac- 
uum expectation values of the two Higgs doublets, mj^ ^ 
are the soft squark masses. Ad are the soft SUSY breaking 
parameters describing the strength of the trilinear scalar 
interactions, and fi is the supersymmetric Higgs (ino) 
mass parameter. Once squark mixing is included, the 
current eigenstates dL,df( become superpositions of the 
mass eigenstates di with the masses toj and the corre- 
sponding mixing angle 9'^ defined as 



TO J 



sin26i'' 



(5) 
(6) 



where to|^, m|^, m\j^ denote the (1, 2), (1, 1), (2, 2) en- 
tries of the mass matrix ^ . 

The squark-neutrino mechanism is a long-range mech- 
anism involving both heavy squark and light neutrino vir- 
tual states as shown in Fig. 1(a). The bottom part of the 
diagram corresponds to the ordinary SM charged current 
interaction. The Ql ~ Qr - mixing results in lepton num- 
ber violation with AL = 2 through the l^p interactions 
proportional to the A'A'-couphng. Without the ql — Qr 
- mixing the l^p -contribution would be proportional to 
the A' A'* -coupling, conserving lepton number. Thus, in 
the latter case LNV with AL = 2 is introduced by the 
Majorana mass of the neutrino with the amplitude of 
Oj^/3/3-decay proportional to the small neutrino mass and, 
therefore, is strongly suppressed. 

It is straightforward to derive the effective 4-fermion 
V — u — d ~ e vertex induced by the squark exchange cor- 
responding to the top part of the diagram in Fig. 1(a). 
The corresponding effective Lagrangian, after a Fierz re- 
arrangement, takes the form 



G 



^{q)LR 



8V2 

[(P, a'^''(l + 75)e^)(wa^.(l-H75)d) 
-t- 2 (p,(l-H75)e^)(u(H-75)d)], (7) 

with (7^'^ = (i/2)[7^, 7*^]. Here U is the neutrino mixing 
matrix and the SUSY LNV parameter is defined as 



nj 
^{q)LR 



y^^i^sin2<, 



1 



1 



2 



d2{k) , 



(8) 



Here we use the notations = d, s, h. This LNV pa- 
rameter vanishes in the absence of ~ Qr. - mixing when 
6'^ = 0, in accordance with the previous comments. 



III. HADRONIZATION PRESCRIPTIONS 

To evaluate the contribution of the diagram in Fig. 
1(a) to Oi^/3/3-decay one first should express the quark 
fields in terms of nucleon ones. This procedure, know 
as hadronization, has so far a rather poor theoretical 
background. In practice the transformation is carried 
out in a phenomenological way, where different possibili- 
ties of embedding the quark fields into the interpolating 
hadronic fields are considered. _ 

There are the two possibilities [9] of hadronization for 
the diagram of Fig. 1 (a) . 

1. The four quark fields can be embedded in the two 
initial neutrons and two final protons separately @, 0j ■ 
This is the conventional two-nucleon (2N) mode of 0i//3/3- 
decay shown in Fig. 2(b). 

2. Another possibility is to incorporate one pair of 
(anti-)quarks of the underlyingquark-level ^pSUSY tran- 
sition into one virtual pion [lol . [Tl| . while another 
(anti-)quark pair participates in the charged current SM 
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(a) 



(b) 



(c) 



FIG. 1: Diagrams describing the squark-neutrino mechanism of Oz^/3/3-decay (a) at the quark-level; at the nucleon level with 
hadronization in the (b) 2N-mode and (c) the pion-mode. The dark blobs denote the LNV vertices generated by I^p SUSY 
trilinear interactions. 



transition involving initial neutron and final proton. In 
this case, which we call the pion mode of Oi^/?/3-decay, the 
second (the upper one in Fig. 1(c)) nucleon vertex is con- 
nected to the leptonic LNV vertex through the charged 
pion exchange. Note that the pion exchange contribution 
in the charged current SM interaction (the bottom vertex 
in Fig. 1(c)) is automatically taken into account by the 
induced pseudoscalar form factor of the nucleon. 

The effective hadronic Lagrangian describing the ver- 
tices of the diagrams in Figs. l(b,c) can be derived using 
the method of Ref. @, matching the hadronic matrix 
elements of the hadronic and the corresponding quark 
operators. This Lagrangian, taking into account both 
the nucleon (p, n) and 7r-meson degrees of freedom in a 
nucleus, has the following form: 



The Lagrangian term Ca in Eq. p4p stands for 
the standard pion-nucleon interaction with the coupling 
g-irNN = 13.4 ± 1 known from experiment. 

The nucleon currents are defined as [l^l 

Jf = <P(p)|u7^(l-75)d|iV(p') > 



P 



2m. 



Jps — 



Jprp 



with 



rPT 



L 



(9) 



(10) 



- 9A{q^)i^ib - 9p{q^)q^ih] n, 

<P{p)\u{l + -i^)d\N{p') > 

FP{q^)+F^;\q^) 
<P{p)\u<j^''{l + ^5)d\N{p') > 



P 



(16) 
(17) 

(18) 



where nip is the nucleon mass, 



rPS 



1^7 



{p — p')^ is the mo- 
mentum transferred to the nucleon vertex with p and p' 
being the four momenta of neutron and proton, respec- 
tively. The tensor structure is given by 



zTf (g2) 



mp 



(19) 



= ^vI^^lrJ'pt E^e>".<^..(l + 75)e^ (12) 



Yqn 



-^vlq)LRmlKn^ Ul,v,{l -I- 75)6^^(13) The nucleon form factors gv{q^), guiq^), gA{q^), gp{q^) 



gT:NN P il^n vr^ 



(14) 



Here, Cp is the charged current SM term, , 
and £^ are the I^p SUSY induced nucleon-electron- 
neutrino and pion-electron-neutrino interactions, respec- 
tively. The coefficient a^r in Eq. (|13p is determined by 
the matrix element Q 



Fg , Fp ,Tj^ are real functions of the squared momen- 
tum q^ transferred to a nucleon. For all the form factors 
we use the following dipole parameterizations: 

gv,AMq^) Fs%(q') Ti'\q') 



gv,A,M 



i|M75d|7r ) = iV2fj, 



rriu + md 



im^K, (15) 



with Ay — 0.71 [GeVY and the normalization constants: 
gv ^l,gA = 1.254, = {np-Hn)gv, {^J■p-^'■n) 3.70. 
The induced pseudoscalar coupling is given by the PCAC 
relation 



where /^r = 0.668 m-^. 



gp{q ^) = 2mpgA{q ^)/{ml ~ q '^). 



(21) 
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TABLE I: Normalizations of nucleon form factors at (f = 
calculated in the quark bag model (QBM) and the non- 
relativistic quark model (NRQM). The table is taken from 
Ref. [Hi. 



Set fP 


y(3) 


rp{2.) 


rp{2.) 


QBM 0.48 4.41 


1.38 


-3.30 


-0.62 


NRQM 0.62 4.65 


1.45 


-1.48 


-0.66 



Using the method of second quantization in the relative 
coordinates we obtain 

=< i/lp_GT(^12)^12 + H%_^{n2)Sl2) >, (24) 

M'' -< 

"^MT ^MT-GT 



{ri2)<Ji2+Hlj,_j,{ri2)Si2) >, (25) 



(26) 



B) The pion mode nuclear matrix element obtained 
within the same framework is 



Af| Hl^_^^{n2)<J,2 + H%_^{n2)S,2) > ■ (27) 



For the normalization constants F^PpiO), T^'^'iO) we use 
the results of Ref. summarized in Table HI obtained 
within the quark bag model and the non-relativistic 
quark model. 



IV. Oi^Pf] HALF-LIFE AND NUCLEAR MATRIX 
ELEMENTS 



The matrix element of the SQN-mechanism (SQuark- 
Neutrino-mechanism) can be calculated according to the 
diagrams of Fig. l(b,c) with the vertices described by 
the effective Lagrangian of Eq. ([9]). The vertices de- 
noted by black blobs originate in squark exchange and 
correspond to the terms pT|) - P^ . The bottom parts 
of these diagrams are the standard model charged cur- 
rent interactions (jlOp . Applying the standard procedure 
based on the non-relativistic impulse approximation (for 
details see [l^), it is straightforward to derive the nuclear 
matrix elements and the corresponding half-life formula 
for the Oi^/3/3-decay in the O"*" 0+ transition channel 
with two outgoing electrons in S-wave states. It takes 
the form: 



Tl/2 



= Goi\Mlf\vl^^^n\ 



(22) 



where Gqi is a precisely calculable phase-space factor. 
Note that it is equal to the phase-space factor of the clas- 
sical light Majorana neutrino exchange mechanism. The 
nuclear matrix elements depend on the hadroniza- 
tion mode. We denote the matrix elements for the 2N 
and pion modes as and A/|, respectively. 

We derive these nuclear matrix elements up to the or- 
der of l/rup in the non-relativistic expansion. They have 
the following form. 

A) The two nucleon mode results in the full matrix 
element 



(23) 



The partial nuclear matrix elements Mj^p, M^^rp and 
Myp have their origin in the interference of the axial- 
vector and pseudoscalar currents, weak-magnetic and 
tensor currents, and vector and tensor currents, respec- 
tively. We note that the Myp contribution was neglected 
in previous studies p^. 



ri2 



ri2 



|ri2|, ri2 



In the above formulas we use the notations: 

£12 

ri2 ' 

CU = (?1 • f?2 

Si2 = 3(cti • fi2)(CT2 • fl2) - 0'12. 



(28) 



Here, ri and r2 are the coordinates of the nucleons un- 
dergoing beta decay. The form of the matrix element 
< 0(1,2) > within the pn-QRPA will be presented in 
the next section. 

The neutrino-exchange potentials in Eqs. ([24)) - (p6)l . 
(P7)) are defined as 



type- 



2 



fK{qri2) - , . r , 



(29) 



where type - K stands for VT - F, MT - GT, MT - 
T, AP - GT, AP -T,nN ~ GT, ttN - T. We also de- 
fine fF..GT{qr 12) = ioiq'T 12) and friqr 12) = -J2(9?'i2), 
where jo,2 are the spherical Bessel functions. In the de- 
nominator of Eq. (|29p E'^g ^ and E^ ^ are the ground 
state energies of the initial and final nuclei, respectively, 
while is the energy of the intermediate nuclear states. 
R = r^A^/^ is the mean nuclear radius with rg = 1.1 fm. 
The functions h1yp^_j^ in Eq. ([29]) can be expressed in 
the following form: 



gv{r){T^'\r)^2Tr>ir) 



(3), 



-• 2 

q 



impTrif, ' 



(30) 



4/T-GT('r') = -3(5y(g') + .9A./(<f^))rf)(g2) 



ivTLpme ' 



/,9 /;?2^ 



"AP-GT 



gAq')Fl.'\q') 



h1pp_rp{q'^) 



2mpme q"^ + ' 
{q% 



(31) 



(32) 



(33) 
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TABLE II; Nuclear matrix elements (NMEs) of the squark-neutrino SUSY mechanism of OufJ/J-decay. The NMEs of the 
2N-mode are calculated for the two cases of the nucleon form factors: Quark Bag Model (QBM) and Non-Relativistic Quark 
Model (NRQM). The quantities M2N, are the 2N and pion mode nuclear matrix elements averaged over small, medium 
and large model spaces (see the text) with their variance a given in parentheses. 







QBM 








NRQM 






nucl. My J, 


M'mt 


MIp 


ML 






MIp 


ML 


M$ 


'^^Ge 185. 


-246. 


29.6 


-20.0 (33.8) 


102. 


-258. 


31.2 - 


■113. (25.7) 


604. (74) 


^°°Mo 220. 


-244. 


33.0 


22.9 (1.8) 


121. 


-256. 


34.8 - 


■100. (14.3) 


594. (80) 


"°Te 179. 


-206. 


28.4 


6.3. (20.8) 


99.2 


-217. 


29.8 ■ 


-81.6 (17.2) 


517. (32) 



hi, 



"ttN-GT 



1 2/-2n^^'4^_ 

R gA[Q — 7 : T 

6 7ne(rnu + md) 



(34) 
(35) 



where the nucleon form factors gK {K = V, A and M), 
F^^") and T'-^'> are defined in sect. III. In the case of 
the pion-exchange mechanism, Eq. we used the 

Goldberger-Treimann relation [13 • 

The following comment is in order. Within the con- 
sidered non-relativistic approximation the NMEs of the 
2N-mode, Eqs (|30 l) -([33 l) . depend nearly on all the nucleon 

(3) (3) 

form factors except for Fg and Tg . As seen from Table 
U] these form factors are rather model dependent quan- 
tities. On the other hand the pion-mode NMEs, Eqs. 
([M)l - (|35p . involve only the axial- vector form factor gj^, 
which is well-known from experimental measurements. 
Thus, the pion-mode NMEs are expected to be signifi- 
cantly less dependent on nucleon structure models than 
the NMEs of the 2N-mode. 



V. CALCULATION OF NUCLEAR MATRIX 
ELEMENTS 

We calculate the NMEs introduced in the previous sec- 
tion. The partial matrix elements in Eqs. ([M)l - ([?71) are 
given by the expression [19|] : 



K 



E E 

J'^ ,ki,kf ,3' pnp'n' 



(36) 



jp jn J 
jn' jp' yJ 



V2J+1 

(p(l),p^2); J II /(ri2)0/^/(ri2) || 7^(1), n'(2); ^) x 
(0+||[c+c„.]j||J-fc/)(J'^fc/|J'^fc.)(J^A;;z||[c+a„],;||0+) . 

The sum over J'" represents the summation over the 
states of the intermediate odd-odd nucleus with angular 
momentum J, parity tt and energies E^jX^' . The indices 
ki and fc/ represent states which are calculated from the 



ground state of the initial (i) and final (/) nuclei, respec- 
tively. The overlap factor {J'^ k f\J'^ ki) accounts for the 
difference between them. The reduced matrix elements 
of the one-body operators c„ (c„ denotes the time- 
reversed state) in depend on the BCS occupation 
coefficients Ui^Vj and on the QRPA vectors X^Y [19j . 

The operators OK,with K = F, GT, T, in ^ depend 
on the distance ri2 between two inicial neutrons that are 
transferred into two protons, the relevant spin and isospin 
operators as well as on the energies of the excited states 
E*yx^' , albeit in practice rather weakly. Short range cor- 
relations of the two initial neutrons and the two final 
protons are described by the Jastrow function f{ri2) ac- 
cording to Ref. pll] . 

We apply the pn-QRPA to calculate the Oi^/3/3-decay 
NMEs for ^^Ge, ^^°Mo and '^^^Te. For these nuclei 
strong experimental limits on the Oj^/3/?-decay half-life 
were found [see Eqs. H]) and ^] and these nuclei are 
also considered as candidate sources for the next gener- 
ation of the 0i^/?/3-decay experiments. For each of them 
three choices of single particle basis (minimal, intermedi- 
ate and large s.p. model spaces) are considered according 
to Ref. [2^. The single particle energies are obtained 
from a spherical. Coulomb corrected Woods-Saxon po- 
tential. The Brueckner G-matrices of the Bonn potential 
are used as a two-body interaction. The pairing inter- 
actions of the nuclear Hamiltonian H are adjusted to fit 
the empirical pairing gaps. In addition, we renormalize 
the particle-particle and particle-hole channels of H by 
introducing the parameters gpp and gph, respectively. In 
the calculation we use gph = 1.0 and gpp is fixed so that 
the known half-life of the 2vB§-(iecsy is correctly repro- 
duced for each model space [23|. Recently, it was found 
that such a procedure makes the Qvl3(3-dec&y NMEs es- 
sentially independent of the size of the single particle 
basis and the nuclear structure input [12, H^l ■ This pro- 
cedure allowed us also to evaluate the uncertainties in 
the calculated NMEs. 

Numerical results of the NME calculation are summa- 
rized in Table [III The presented values of the partial nu- 
clear matrix elements Myj,, 



Mljrp and M\p have been 
obtained for the intermediate s.p. model space: 12 levels 
for A=76, 16 levels for A=100, and 18 levels for A=130. 
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The values of NMEs for the three studied isotopes are 
comparable in magnitude. We have found that the con- 
sidered two-nucleon short range correlations suppress the 
pion-mode NME by about 20%. We recall that the sim- 
ilar reduction occurs also for the NMEs of the "classic" 
light Majorana neutrino mass mechanism of the 
decay. 

As seen from Table |lTl the pion-mode NMEs are from 
6 to 30 times larger, depending on the nucleon structure 
model, than the corresponding 2N-mode NMEs for all 
the studied nuclei. This allows us to conclude that in the 
squark-neutrino I^p SUSY mechanism of Oz//3/3-decay the 
pion-mode of hadronization clearly dominates over the 
2N-mode. 



TABLE III: Upper bounds on the |tpSUSY parameter r)\l^^ji 
as well as on the related products of the trilinear _^p-couplings 
-^'iifc-^'ifci (k=li2,3) for KsusY = 100 GeV (see scaling law in 
Eq. ((STjl) deduced from the current lower bounds on the half- 
life of Oi^/3/3-decay for "^^Ge, ^°°Mo and ""Te. 



nucl. 


T";--'' [Ref.] vr,)Lu 
(years) 


-'^iii-'^iii 


-^112-^121 


'^113-^131 




> 1.9 10^'' [2| 4.3 10"'^ 


7.7 10"'' 


4.0 10"'' 


1.7 10"* 


Mo 


> 5.8 10^^ [4] 9.2 10"'' 


1.7 10"^ 


8.7 lO"'' 


3.6 10"® 




> 3.0 10^* f5] 4.7 10"^ 


8.5 10"'' 


4.5 10"^ 


1.9 10~® 



Let us point out that in the 27V-mode NMEs there is 
a strong cancelation between the Myj, and M^^jj, partial 
matrix elements (see Table |TT]) originating from the tensor 
nucleon current As a result despite the individual 

values of both matrix elements Myj, and M^jj, are sig- 
nificantly larger than the pseudoscalar contribution M\p 
it turns out that their cumulative effect, Myj, -\- Mf^j,, is 
suppressed and can be either comparable in magnitude 
with Af^p or larger than it, depending on the model of 
nucleon structure (QBM or NRQM). At this point we dis- 
agree with Ref. tl5(] claiming that the tensor contribution 
is always significantly larger than the pseudoscalar one. 
This disagreement can be accounted for the fact that the 
contribution M^j, was missed in Ref. [l5| leaving the 

large term, analogous to our M^jj,, uncompensated. As 
a result the tensor contribution, represented by only this 
term, becomes significantly larger than the pseudoscalar 
contribution for both QBM and NRQM models of nu- 
cleon form factors. As we have shown this is not the case 
when all the contributions of the tensor nucleon current 
(1181) are properly taken into account. 

In Table [U the NMEs of the pion and 2N modes are 
give n with their variances a (evaluated according to Ref. 
|23 |) due to uncertainties of nuclear structure input. In 
both cases these variances are within 10%. This compar- 
atively low uncertainty was obtained mainly due to the 
procedure we chose for fixing the residual interaction of 
the nuclear Hamiltonian which includes the calculation 
with three model spaces [13, HH . As to the uncertainties 
of the nucleon form factors, from Table [ll] it is evident 
that they are very significant for the 2N-mode NMEs. For 
the two models of nucleon form factors considered here, 
the Quark Bag Model and the Non-Relativistic Quark 
Model, the values of the partial NMEs of the 2N-mode 
change up to a factor 2 while the total 2N-mode NMEs 
up to an order of magnitude. Fortunately, as we com- 
mented at the end of Sect. IV the dominant pion-mode 
contribution is free of the nucleon structure uncertainties 
making their impact on the resulting NMEs insignificant. 



VI. DISCUSSION AND CONCLUSIONS 

With the values of the nuclear matrix elements of the 
dominant pion-mode contribution given in Table [Til we 
derived the experimental upper limits on the LNV pa- 
rameter r]lq)LR introduced in Eq. ([7]). Then, using Eq. 
([5]), we extracted upper limits on the products of the 
trilinear l^p -couplings ^n^^iki (k=l,2,3). These lim- 
its, presented in Table IIIIl have been derived under the 
conventional simplifying assumptions. We assumed all 
the squark masses and the trilinear soft SUSY breaking 
parameters Ad to be approximately equal to a common 
SUSY breaking scale Asusy- Thus we approximately 
write 




with gfe = (1.8 X 10^; 94.2; 3.9) calculated for the cur- 
rent quark masses nid — 9 MeV, nis = 175 MeV and 
nib = 4.2 GeV respectively. In the above formula T°^2~^^^ 
is the experimental lower bound on the half-life of a cer- 
tain nucleus. In Table lllll we show the upper limits on the 
SUSY parameters, extracted from the existing exper- 
imental lower bounds on the half-lives of the three nuclei 
™Ge, i°°Mo and ^^°Te. Our limit for '^^Ge is about an 
order of magnitude better than the limit previously ob- 
tained in Ref. [ij, [lH| on the basis of the 2N-mode of 
hadronization taking into account only the pseudoscalar 
current contribution, M'^p. 

In conclusion, we analyzed Oi^/3/3-decay of several nu- 
clei induced by the LNV effective operators originating 
from the I^p SUSY trilinear interactions involving squark 
and neutrino exchange. We focussed on the hadroniza- 
tion prescription of the quark-level operators and ana- 
lyzed both the conventional 2N-mode and the pion-mode 
of hadronization. We have shown that the pion-mode 
absolutely dominates over the 2N-mode. Previously, in 
Refs. [1] it was demonstrated that the pion-mode dom- 
inates over the 2N-mode in the case of the short-range 
-^pSUSY mechanism. Thus, with the result of the present 
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paper we conclude that all the mechanisms based on the 
trilinear SUSY interactions dominantly contribute to 
0://3/3-decay via the pion-mode of hadronization. 
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